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Summary 

The vertical wind or downdraft component of 
a microburst- generated wind shear can significantly 
degrade airplane performance. Doppler radar and 
lidar are two sensor technologies being tested to pro- 
vide flight crews with early warning of the presence 
of hazardous wind shear. An inherent limitation of 
Doppler-based sensors is the inability to measure ve- 
locities perpendicular to the line of sight; this results 
in an underestimate of the total wind shear hazard. 
One solution to the line-of-sight limitation is to use 
a vertical wind model to estimate the vertical com- 
ponent from the horizontal wind measurement. The 
objective of this study was to assess the ability of 
simple vertical wind models to improve the hazard 
prediction capability of an airborne Doppler sensor 
in a realistic microburst environment. Both simula- 
tion and flight test measurements were used to test 
the vertical wind models. The results indicate that, 
in the altitude region of interest (at or below 300 m), 
the simple vertical wind models improved the hazard 
estimate. The radar simulation study showed that 
the magnitude of the performance improvement was 
altitude dependent. The altitude of maximum per- 
formance improvement occurred at about 300 m. At 
lower altitudes, the improvement was minimized by 
the diminished contribution of the vertical wind. The 
vertical hazard estimate errors from flight tests were 
less than those of the radar simulation study. 

Introduction 

Wind shear is considered by many in the avia- 
tion industry to be a major safety issue. Numerous 
accidents and incidents have occurred that were at- 
tributed to low-altitude wind shear, which can be 
found in a variety of weather conditions such as gust 
fronts, sea breeze fronts, and mountain waves (ref. 1). 
However, hazardous wind shear is most often associ- 
ated with the convective outflow of thunderstorms 
known as microbursts. A microburst is a strong lo- 
calized downdraft that causes a significant outflow 
as it impacts the ground (ref. 2). The hazard of 
a microburst encounter occurs when a head wind 
rapidly shifts to a tail wind as the airplane pene- 
trates the outflow, which reduces the airspeed and 
the potential rate of climb of the airplane. The po- 
tential rate of climb of the airplane is further reduced 
by the microburst downdraft. The general effect on 
the airplane is a rapid loss of energy from which it 
may not have enough altitude, airspeed, or thrust to 
overcome. 

The National Aeronautics and Space Adminis- 
tration, in a joint effort with the Federal Aviation 
Administration, has been conducting research in the 


development of forward-look, airborne, wind shear 
detection technologies. Forward-look systems warn 
the flight crew of the presence of wind shear at low 
altitude (under 305 m) in time to avoid the affected 
area or to prepare for and escape from the encounter. 
A fundamental requirement of such a wind shear de- 
tection system is the ability to reliably estimate the 
magnitude of an upcoming wind shear hazard along 
the flight path of the airplane. Doppler radar and li- 
dar are two technologies being tested to provide this 
capability. Both measure the Doppler shift from rain 
drops, aerosols, and other debris in the air to deter- 
mine the line-of-sight relative velocity of the air. 

An inherent limitation of this system is its inabil- 
ity to measure velocities perpendicular to the line 
of sight. The presence of a microburst can be de- 
tected by measuring the divergence of the horizontal 
velocity profile; however, the inability to measure the 
downdraft can result in a significant underestimate 
of the magnitude and spatial extent of the hazard 
(ref. 3). 

One solution to the line-of-sight limitation of 
Doppler sensors is to use a theoretical or empirical 
model of a microburst to estimate the perpendicular 
velocities from the measured line-of-sight values. A 
preliminary assessment of this technique showed that 
the microburst-generated downdraft could be esti- 
mated at low altitudes with very simple vertical wind 
models (ref. 3). This study was, however, limited in 
scope. It assumed perfect knowledge of the horizon- 
tal wind field produced by an idealized axisymmetric 
microburst simulation. 

The objective of this study was to assess the 
ability of simple vertical wind models to improve the 
hazard prediction of an airborne Doppler sensor in 
a realistic microburst environment. Both simulation 
and flight test measurements were used to test the 
vertical wind models. A Doppler radar simulation 
was used with a high-fidelity asymmetric microburst 
model to establish the performance limits of the 
vertical wind models and to establish the effects of 
radar signal noise and measurement errors. Flight 
test measurements from an airborne Doppler radar 
were used to compute the vertical wind velocity in 
front of the airplane. These velocities were compared 
with onboard in situ measurements. The flight test 
results were also compared with the expected results 
established from the radar simulation. 

Although the simulation and the flight test re- 
sults presented in this report are for an airborne 
Doppler radar system, the techniques for estimat- 
ing the vertical wind velocity should be applicable to 
other Doppler-based sensors such as lidar. However, 



the performance will vary with the signal properties 
of the sensor and the spatial resolution. 

This paper reviews the wind shear hazard index 
known as the F- factor, which is used extensively to 
quantify the hazard of a wind shear encounter. The 
microburst simulation data sets used in the radar 
simulation are introduced. The coordinate system, 
the two vertical wind velocity estimation techniques, 
the radar simulation, and the flight test data are then 
described. The method of analysis and results are 
also presented. 

Symbols 

F wind shear hazard index 

Ffr horizontal component of wind shear 

hazard index 

F v vertical component of wind shear 

hazard index 

F wind shear hazard index averaged over 

specified distance 

F v vertical wind shear hazard index 

averaged over specified distance 

/(r m ) empirical model radial shaping 

function of radial wind velocity, m/s 

g gravitational acceleration, m/s 2 

#( r m) empirical model radial shaping 

function of vertical wind velocity 

p{z m ) empirical model vertical shaping 

function of radial wind velocity 

q(z m ) empirical model vertical shaping 

function of vertical wind velocity, m/s 

R linear correlation coefficient 

Res residual of linear least-squares curve 

fit of radial velocity profile, m/s 

r c sensor radial coordinate to origin of 

microburst reference system, m 

r m microburst-referenced radial 

coordinate, m 

r ma x microburst radial coordinate of 

maximum horizontal wind, m 

r min minimum radar range 

r s airplaine sensor- referenced radial 

coordinate, m 

T^ 5 velocity transformation matrix for 

sensor elevation angle 


Tfi m velocity transformation matrix for 

microburst azimuth angle 

T# s velocity transformation matrix for 

sensor azimuth angle 

t time 

U m microburst-referenced wind 

vector, m/s 

U s sensor-referenced wind vector, m/s 

Uqo microburst- referenced free-stream wind 

vector, m/s 

u m microburst- referenced radial wind, m/s 

u s sensor- referenced radial wind, m/s 

Uqo microburst-referenced free-stream 

radial wind, m/s 

it time rate of change of horizontal wind 

component (tailwind positive), m/s 2 

V true airspeed at time of radar 

measurement, m/s 

V g ground speed, m/s 

v m microburst-referenced velocity in 0 m 

direction, m/s 

v s sensor-referenced velocity in 8 S 

direction, m/s 

w vertical wind component (updraft 

positive), m/s 

w s sensor- referenced velocity in 4> s 

direction, m/s 

w m microburst-referenced vertical wind 

component, m/s 

x coordinate in west-east direction about 

center of microburst, m 

y coordinate in south-north direction 

about center of microburst, m 

z coordinate in vertical direction, m 

z m microburst-referenced vertical 

coordinate, m 

Zmax altitude of maximum horizontal 

wind, m 

z s altitude of sensor, m 

a empirical model shaping function 

variable 

A empirical model scaling factor, s -1 
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A r s sensor range bin length, m 

At time shift between radar and in situ 

measurement, s 

7} (z m ) empirical model altitude weighting 

function, m 

0 C sensor azimuth coordinate to 

microburst reference system 
origin, deg 

6 m microburst-referenced azimuth 

coordinate, deg 

9 S airplane sensor-referenced azimuth 

angle, deg 

0 o o microburst-referenced azimuth of free- 

stream wind, deg 

4>s airplane sensor- referenced elevation 

angle, deg 


The F- factor can be separated into a horizontal 
component F ^ and a vertical component 1 F v , so that 


where 


? = Fh + Fv 

(2) 

F “ 

*h i — ~ 

(3) 

9 

w 

Fy = 

V 

(4) 


Doppler-based wind shear sensors can only mea- 
sure the line-of-sight divergence of the wind and 
therefore can only determine the horizontal F-factor 
component. The inability to measure the vertical 
component can result in a significant underestimate 
of the magnitude of the microburst hazard. For land- 
ing or departing airplanes, the vertical component of 
F can exceed half the total F-factor. 


Microburst Simulation Model 


Abbreviations: 


AWDRS 

TASS 

2D 

3D 


Airborne Windshear Doppler Radar 
Simulation Program 

Terminal Area Simulation System 

two-dimensional 

three-dimensional 


Wind Shear Hazard Index 

The magnitude of the hazard posed by a 
microburst to an airplane can be quantified through 
the F-factor (ref. 4). The F-factor is a hazard in- 
dex that represents the rate of specific energy loss 
because of wind shear. For straight and level flight, 
the F-factor can be expressed as 


Positive values of F indicate a performance- 
decreasing condition. Conversely, negative values 
indicate a performance-increasing condition. The 
F-factor is directly related to the climb gradient ca- 
pability of the airplane. For example, a value of F of 
0.2 would indicate a loss in climb gradient capability 
of 0.2 rad (11.5°). If an airplane had a maximum 
climb angle capability of 10° , it would be unable to 
maintain level flight in that sustained wind shear en- 
vironment. A wind shear is considered hazardous to 
landing or departing airplanes if the F-factor average 
over 1 km is 0.1 or greater (ref. 4). 


The Terminal Area Simulation System (TASS) 
high-fidelity microburst simulation model was used 
in the development and the analysis of the down- 
draft estimation techniques discussed in this paper. 
TASS is a time-dependent, multidimensional, non- 
hydrostatic, numerical cloud model that has been 
used extensively in the study of microbursts (refs. 5 
and 6). The model is initiated with the ob- 
served environmental conditions (altitude profiles 
of temperature, humidity, and wind) that existed 
prior to microburst development and outputs a 
three-dimensional time history of radar reflectivity, 
winds, temperature, pressure, water vapor, rainwa- 
ter, snow, hail, and cloud water. The model has 
been validated with both ground-based and airborne 
measurements (refs. 6 and 7). Microburst data sets 
generated with TASS have also been selected by the 
Federal Aviation Administration as test cases for cer- 
tifying forward-look wind shear sensors (ref. 8). Al- 
though TASS is much too complex to be practical 
as a downdraft estimation model, it is very useful 
for generating the high-fidelity data sets necessary 
to evaluate such models. 

Two TASS microburst data sets were used in this 
study. One was an axisymmetric case (symmetric 
about the vertical axis) with a 20-m grid resolution, 
and the other was a three-dimensional asymmetric 
case with a 200-m horizontal grid. 

Axisymmetric Microburst Data Set 

The axisymmetric data set used in the study 
reported herein was also used in an earlier feasibility 
study of vertical wind estimation techniques (ref. 3). 
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In this study, it is used to illustrate the general 
characteristics of a microburst wind field. These 
characteristics are discussed in the section “Vertical 
Wind Models.” The axisymmetric data set extended 
from the microburst core to 4000 m radially and 
from the ground to 600 m vertically. This was a 
subset of the larger TASS modeled domain. The 
model was initiated with the atmospheric conditions 
measured before a microburst event on August 2, 
1985, at Dallas-Fort Worth International Airport 
(ref. 9). Four different times in the microburst 
simulation were selected, each 2 min apart. Figure 1 
shows the wind vector plots for the four time periods 
selected. The first was just before the downdraft 
impacted the ground, which was at 9 min into the 
microburst simulation. The second was just after 
the downdraft hit the ground and began to spread 
out, which was approximately the time of maximum 
horizontal shear. The third was at a point when the 
outflow vortex was well defined, and the last, near 
the end of the life cycle of the microburst event. 

Contour plots of F-factor for an airplane flying 
level at 130 knots are also shown in figure 1. Figure 2 
shows the same data sets with the F-factor contours 
computed without the vertical winds. The magni- 
tude and spatial extent of the detectable hazard are 
clearly diminished. This further illustrates the need 
for some means of determining the magnitude of the 
vertical winds. 

Asymmetric Microburst Data Set 

The asymmetric data set was used in the Doppler 
radar simulation, which is described in section 
“Radar Simulation.” The data set was gener- 
ated from an atmospheric sounding taken before a 
rnicroburst event on July 11, 1988, in Denver. This 
event was inadvertently encountered by four suc- 
cessive airplanes on final approach to the Stapleton 
International Airport in Denver. This data set accu- 
rately simulates the major features of the microburst- 
producing storm and is discussed in detail in 
reference 7. The data set extends 12 km in the 
west-east direction (^-direction), 12 km in the south- 
north direction (^-direction), and 2 km vertically 
( z-direction), with a vertical resolution of approx- 
imately 80 m. Only the data set from the lowest 
600 m was used in this study. Three times from this 
simulation were selected, each 1 min apart. Figure 3 
shows the horizontal cross section of the wind vector 
field at an altitude of 283 m for each of the three 
times selected. Also shown in the figure is an out- 
line of the simulated radar scan area that is discussed 
later. At the bottom of the figure is the vertical cross 
section of the wind vector field through the location 


y — 0. The three times selected correspond to an in- 
stant roughly 1 min prior, during, and 1 min after the 
first of the four airliners encountered the microburst. 
These times are also near the time of maximum shear 
produce by the microburst. 

Velocity Transformation Equations 

Estimating the vertical wind with a microburst 
model requires a transformation between the 
cylindrical coordinate system (r m ,0 m ,2 m ) of the 
microburst and the spherical coordinate system 
( rs^sfis ) of the airplane sensor. The sensor is as- 
sumed to be inertially referenced so that the airplane 
motion is compensated for and is removed from the 
sensor- measured velocities. Figure 4 shows the two 
coordinate systems and their corresponding veloc- 
ity components. The wind field is assumed to be 
a microburst wind field superimposed on a uniform 
horizontal wind field. The microburst vertical axis is 
located in the sensor coordinate system at (r c ,0s,0 c ). 
The origin of the sensor coordinate system is at an 
altitude z s above the origin of the microburst co- 
ordinate system. Derivation of the velocity trans- 
formation equations and their spatial derivatives are 
provided in the appendix. 

Vertical Wind Models 

Two vertical wind models are analyzed in this 
report. The derivation of these models is provided 
in this section. Both of the models determine the 
vertical divergence dw m /dz m from the measured ra- 
dial wind profile with the mass continuity equation. 
The vertical wind w m is then determined from a 
model-based relationship between the vertical diver- 
gence and the vertical wind with altitude. 

Estimating Vertical Divergence 

The mass continuity equation for incompressible 
flow in the microburst cylindrical coordinate system 
is 

1 d(r m u m ) 1 dv m dw rri /-v 

7 1 ™ r- — u V°y 

Tjyi OTjxi T fa UUm uZm 

If the microburst is assumed to be symmetrical about 
the vertical axis (dU m/d9 m = 0) with no rotational 
flow (v m = 0), then equation (5) simplifies to 

Um . dw m /g\ 

r m dr m dz m 

The microburst radial velocity u m and the radial 
shear dum/drm can be determined from the sensor 
measurements and the transformation equations of 
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the appendix. If the sensor elevation angle is assumed 
to be level with the horizon (<f> s = 0), then equa- 
tion (7), which is equation (A46), can be used to 
relate the sensor- measured radial shear du s /dr s to 
the microburst-referenced values: 


dug 

dr s 


dum 

d r m 


cos 2 (& m - 0 S ) + — sin 2 (0 m - 0 S ) (7) 

r m 


Figures 5 and 6, reproduced from reference 3, 
show the radial profiles of the radial and the verti- 
cal wind from the symmetrical microburst data set. 
Figure 5 shows that near the core of the microburst 
(from r = 0 to r ~ 1200 m) the radial velocity varia- 
tion is nearly linear. Comparing figure 5 with figure 6 
shows that this linear region corresponds to the pri- 
mary downdraft region of the microburst for each of 
the four times. The vertical velocity variations out- 
side the core area are due to the outflow vortex ring, 
which expands radially with time. In the linear core 
region, 

dUffi Urn 

^ « — (8) 

or m r m ' 

and equation (7) simplifies to 

dlLm . . 

s Tm 

Combining equation (9) with equation (6) yields a 
simple relationship for the vertical wind divergence 
as a function of the sensor- measured radial shear in 
the core of the microburst: 


dWm 

& z m 


f & u 7n A __ ^ 

\ I'm / dr$ 


(10) 


This relationship is not valid outside the microburst 
core where the radial velocity profile becomes non- 
linear. As the distance from the microburst core in- 
creases, the ratio u m /r m becomes small. If this term 
is neglected, then equations (6) and (7) become 


dum dw m 

9v m dzjYi 


( 11 ) 


and 


dus 

dr s 


m 


cos 2 {o m - o s ) 


( 12 ) 


respectively. Since 0 < cos 2 (0 m - 0 S ) < 1, then 


0u s 

dr s ~ dr m 

du s < dw m 
dr s dz m 


(13) 

(14) 


If this inequality is approximated as an equality, 
then the magnitude of the resultant vertical wind di- 
vergence estimate outside the microburst core may 
be low. Underestimating the magnitude of the di- 
vergence outside the core should not adversely im- 
pact the hazard estimate since the vertical wind 
in this region is generally characterized as either 
performance-increasing updrafts or small scale down- 
drafts, as shown in figure 6. Therefore, a practical 
vertical divergence estimate outside the microburst 
core is 

t)u?jn, du c 

< 15 > 

The vertical wind divergence can be estimated 
from the sensor-measured radial wind profile with 
equations (10) and (15). However, these equations 
require the core of the microburst to be identi- 
fied. The method used to identify the core of the 
microburst from the sensor- measured radial veloc- 
ity profile is discussed in the section “Analysis and 
Results. ” 


Linear Model 


The linear model is the simplest of the two ver- 
tical wind models used in this study. The model is 
based on the assumption that there is no vertical 
wind at the ground and it increases linearly with al- 
titude (i.e., dw m /z rn = Constant). The vertical wind 
can be computed from the vertical wind divergence 


as 


w m 


< 9^771 

Zrn 

uz m 


( 16 ) 


Figure 7 shows the vertical wind variation with 
altitude at the center of the microburst presented in 
figures 5 and 6. The assumption of linearity appears 
reasonable near the center, particularly at altitudes 
below 400 m. At higher altitudes the linearity begins 
to break down. Figure 8 shows the vertical wind 
profile at a radius of 2000 m; the linearity breaks 
down at this location as well. This nonlinearity 
occurs primarily because of the outflow vortex, which 
is generated when the downdraft core impacts the 
ground and begins to diverge horizontally (ref. 3). 
The degree to which these nonlinearities introduce 
errors into the downdraft calculation is discussed in 
the section “Analysis and Results.” 


Empirical Model 

The empirical vertical wind model is based on 
measurements of several microburst events. This 
model utilizes the empirical microburst model devel- 
oped by Oseguera and Bowles (ref. 10) and subse- 
quently modified by Vicroy (ref. 11). The model is 
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an axisymmetric, steady state model that uses shap- 
ing functions to simulate boundary layer effects and 
to satisfy the mass continuity equation. The mass 
continuity equation (eq. 6) can be satisfied by solu- 
tions of the form 


v-m = f(r m )p(Zm) (17) 


The values of c x , c 2 , and z ma x were selected from 
curve fits of data from several microburst events. 
Differentiating equation (18) with respect to altitude 


yields 


dw r 


0z„ 


9iTm) 


2 9q{z m ) 


dzr, 


(28) 


Substituting equation (20) for dq(zm)/dz m yields 


Wm = g{rm)q{zm ) 

provided 

d[r m f(r m )} _ A A 
a 4 - 2 9(rm) 

^ = — Ap(zm) 
dz-m 

where /(r m ) and g(r^) are the empirical radial shap- 
ing functions for the radial and the vertical wind ve- 
locities, respectively; p[z m ) and q(z m ) are the ver- 
tical empirical shaping functions for the radial and 
the vertical wind velocities, respectively; and A is a 
scale factor. The characteristic shape of the radial 
and the vertical shaping functions is shown in fig- 
ures 9 and 10. The shaping functions are used to ap- 
proximate the characteristic profile of the microburst 
winds. The radial shaping functions appear to com- 
pare well with the axisymmetric microburst profiles 
presented in figures 5 and 6, particularly for the first 
two times. The radial profiles of the last two times 
show the growth of the outflow vortex, which is not 
modeled by the shaping function. 

The equations for the shaping functions are 

2 - ( r 2 /r 2 V*' 

“ V mi ' max/ /oil 

2a V 


,, , Ar m 

firm) = exp 


(18) 

(19) 

(20) 


9( r m) = 


1 - 


exp 


9 (v^ f 

L V 771/ ' III ax ^ 


p(zm) — exp ^Ci • 


z m \ 
2 max / 


exp i C 2 


2a 


Zm \ 


Z max / 


( 22 ) 

(23) 


q{ z m) — ^ 


Zmax ( 

exp I Ci - 


Cl 


Zmax 

c 2 


exp C 2 


Zmax 
Zm 


(c2--l) 

V Z ma x / . 


with 


C1 = -0.15 
c 2 = —3.2175 

Zmax ~ hO m 


(24) 


(25) 

(26) 
(27) 


dlVm 

dz m 


— Ap(zm)fl( r m) 


(29) 


Combining equation (29) with equation (18) yields 


-g(-Zm) dwm 
Ap(zm) 9z m 


- T){z m ) 


dWm 

dz m 


(30) 


where 


Zmax 

Cl 


exp ^ci 


Zm 

Zmax 


- 1 


v{Zm) “ ( z m \ ( z m \ 

exp ci 1 - exp I c 2 I 

\ Zmax / V 2 max / 

exp( 


Zmax 


Zm 


- 1 


Zmax 


exp I ci 


( c 


Zrn \ 


Zmax / 


- exp 


(C2 — ) 
\ Z ma x / 


(31) 


The vertical wind can now be computed from the 
vertical wind divergence and the altitude dependent 
function r]. 


Radar Simulation 

The effects of measurement noise and signal pro- 
cessing techniques on the vertical wind estimate were 
determined with the Airborne Windshear Doppler 
Radar Simulation (AWDRS) Program (ref. 12). The 
simulation program computes the airborne radar 
signal returns for a given TASS microburst simu- 
lation and ground clutter map. The radar simula- 
tion includes algorithms for signal filtering and auto- 
matic gain control, as well as computing in-phase and 
quadrature base-band signal components, Doppler 
velocity, spectral width, and F- factor. The AWDRS 
program requires the relative location of the airplane 
and the microburst to the runway touchdown point 
and the radar system parameters, such as antenna 
pattern, pulse width, transmitted power, and an- 
tenna elevation. 

The AWDRS program was initialized with the 
input data shown in table 1. The input data that 
varied between runs, namely the airplane starting 
range and the microburst position, are presented in 
table 2. The TASS generated asymmetric microburst 
data set was used for the radar simulation runs. 
The airplane and the microburst starting positions 
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were selected so that the airplane was at the same 
horizontal location relative to the microburst, and 
at the desired test altitude during the third scan of 
the simulation run. The data from the third scan 
of each simulation were used in the analysis. The 
third scan was selected to allow temporal filtering 
schemes within the simulation to initialize. Scan data 
were collected at six altitudes, from 100 to 600 m 
in 100-m increments, for each of the three times 
(t = 49, 50, and 51 min). The microburst wind field 
was frozen in time during each radar simulation run. 
The radar simulation scan area is outlined on the 
microburst wind vector plots shown in figure 3. The 
azimuth scan covered —21° to 21°, in 3° increments. 
Each scan line consisted of 30 range bins each 150 m 
long, with the initial range bin 425 m in front of 
the airplane. The resultant scan measurement grid 
is depicted in figure 11. The radar antenna elevation 
angle was set to 0 (level with the horizon) for all but 
the two lowest scan altitudes. At those altitudes, 
the elevation angle was set slightly above the horizon 
to minimize ground clutter returns. The elevation 
angle for the 100- and 200- m scans was 1.185° and 

0.470°, respectively. Contour plots of the true radial 
velocity and the simulated radar measurement with 
clutter and noise are shown in figure 12, for each of 
the three microburst cases. 

Flight Test Data 

A series of flight tests were conducted during the 
summers of 1991 and 1992 with a Transport Sys- 
tems Research Vehicle (TSRV) Boeing 737-100 air- 
plane at Langley that was equipped with a vari- 
ety of prototype wind shear detection systems. The 
tests were conducted near Orlando, Florida, and 
Denver, Colorado. These two locations were se- 
lected based on their climatology, which is conducive 
to microburst development, and the availability of 
ground-based Doppler weather radar coverage. The 
Doppler weather radar was used to identify and to 
direct the research airplane to potential microburst 
activity. The microburst penetrations were typically 
flown at airspeeds between 210 and 230 knots and al- 
titudes between 244 and 335 m. Details of the flight 
test procedure can be found in reference 13. 

There were three forward-look wind shear detec- 
tion sensors onboard the airplane: a passive infrared 
sensor, a lidar, and a Doppler radar (refs. 14-17). 
The ground-based Doppler weather radar wind di- 
vergence measurements were also transmitted to the 
airplane via radio data link and processed to compute 
F-factor estimates. The airplane was also equipped 
with a reactive, or in situ, system that computed 
the F-factor of the airspace the airplane was cur- 


rently flying through (ref. 18). The in situ F-factor 
was used as a truth measurement for validation of 
the forward-look wind shear detection sensors. The 
F-factor predicted from the forward- look sensor was 
compared with the in situ measurement of the air- 
plane as it penetrated the scanned airspace. The 
flight data in this study were limited to the onboard 
Doppler radar and the in situ measurements. 

A limited subset of the vast flight data collected 
was selected for analysis in this study. This data 
selection was based on the following criteria: 

1. The Doppler radar had to be operating at 
an antenna elevation of 0 (i.e., level with the 
horizon) with a pulse width of 0.96 (i sec. 

2. The airplane had to be flying nearly straight 
and level. 

These flight requirements ensured that the airplane 
flew through the same airspace that was scanned by 
the radar. The pulse width restriction ensured that 
the radar range bin size was the same bin size used 
in the earlier radar simulation study. The microburst 
events selected for analysis are listed in table 3. 

Analysis and Results 

The analysis of the vertical wind estimation tech- 
niques was conducted in three parts. The objec- 
tive of the first part was to determine how well the 
simple vertical wind models estimate the microburst 
winds independent of signal noise or measurement 
error. This objective established the upper perfor- 
mance limit for the vertical wind models. The second 
part of the analysis focused on the effects of signal 
noise and the measurement error on the vertical wind 
estimation. The TASS generated microburst data 
sets were used as the reference or truth for the first 
and second parts of the analysis. The third part of 
the analysis used the flight test measurements from 
the airborne Doppler radar to compute the vertical 
winds and to compare with the in situ measurements. 
The flight test results were then compared with the 
expected results established in the second part of the 
analysis. 

All the simulation analysis was conducted for alti- 
tudes from 100 to 600 m. This range covers altitudes 
above those for which a forward-look wind shear sys- 
tem would be active. The altitudes above 300 m 
were included in the analysis for comparison with the 
results of the preliminary study of reference 3 and 
to establish the altitude limit of the vertical wind 
models. 
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Validity of Vertical Wind Model 

Assumptions 

The two vertical wind models of this study re- 
late the horizontal shear to the vertical wind as a 
function of altitude. These models share two basic 
characteristics: 

1. Divergent radial winds (positive shear) are 
associated with downdrafts and conversely, 
convergent radial winds (negative shear) are 
associated with updrafts. 

2. The magnitude of the vertical wind is propor- 
tional to the magnitude of the shear. 

An earlier study (ref. 3) established that these 
two assumptions are reasonable for the simple ax- 
isymmetric microburst case presented in figure 1. 
Although some microbursts are roughly axisymmet- 
ric, many are asymmetric. The validity of these 
assumptions has not been established for 

asymmetric microbursts. Figure 13 shows where 

these assumptions are and are not valid for the asym- 
metric microburst data set of this study. Shown in 
the figure is the vertical wind at each radar range 
bin plotted as a function of the radial shear. Also 
shown in the figure are the two vertical wind model 
functions. The data points that lie in the lower left 
and the upper right quadrants of the graphs violate 
the first of the assumed characteristics. The nonzero 
data points that lie along the vertical axis violate 
the second assumed characteristic. Both basic as- 
sumptions appear reasonable at the lower altitudes 
(z < 300) and for the earliest of the three times. 
At the later times and the higher altitudes these as- 
sumptions break down. There appears to be very 
little correlation between the radial shear and the 
vertical wind above 300 m for t — 51 min. 

Microburst Core Criteria 

The vertical wind divergence (which is used in 
the vertical wind models) can be estimated from 
the sensor-measured radial wind with equation (10) 
or (15), depending on whether the estimate is inside 
or outside the downdraft core. Near the downdraft 
core, the radial velocity variation is nearly linear for 
the symmetrical data in figures 5 and 6. Figure 14 
shows the vertical wind contours for the asymmetric 
microburst data at each scan altitude, with the w = 0 
contour line highlighted in white. The white contour 
line was considered to outline the downdraft core. 
On the right side of the figure are the corresponding 
contour plots of the linear correlation coefficient R 
of the radial wind profile. The linear correlation 
coefficient was computed from a moving, 5-point 
linear least-squares fit of the radial wind profile. The 


linear correlation coefficient for range bin i along scan 
line j was defined as 


Rij — 


-2 u*i_2j - + 2u 


s i+2,j 




(32) 


k=i + 2 


io £ u l kJ - 2 1 £ u s , 


k=i — 2 




The value of R varied from —1 to 1 with the magni- 
tude indicating how linear the 5 data points were, a 
value of 1 or —1 being linear and 0 being nonlinear. 
The sign of R corresponded to the sign of the slope 
or shear. 

The w = 0 contour line was repeated on the R 
contour plot to evaluate the correlation between a 
positive linear shear and the downdraft core. As can 
be seen in figure 14, the region of positive linear shear 
encompasses most of the downdraft region. Once 
again, the lower altitudes provided the best results. 
At the higher scan altitudes, there were small regions 
of strong downdraft with either nonlinear or negative 
shears. This condition was assumed to be due to 
smaller microbursts being generated within the larger 
one. 

The effect of measurement noise is shown in fig- 
ure 15, which is the same as figure 14 except that the 
linear correlation coefficient was computed from the 
simulated radar- measured velocities. The measure- 
ment noise and the signal filtering reduce the region 
of positive linear shear, particularly at the lower al- 
titudes where ground clutter is a major factor. 

The correlation between the downdraft core area 
and the area of positive linear shear was signifi- 
cantly eroded when the effects of noise were in- 
cluded. Lacking any other mechanism to distinguish 
the downdraft core, R > 0.9 was selected to represent 
the microburst core area and R < 0.9 the area out- 
side microburst core. 


Radar Simulation Results 

The vertical wind models were first tested inde- 
pendent of signal noise or measurement error to es- 
tablish their upper performance limit. The true ra- 
dial velocities, shown on the left side of figure 12, 
were used to compute the radial shear. The shear 
was computed in conjunction with the linear correla- 
tion coefficient by using a 5-point linear least-squares 
fit. The equation for the radial shear (slope of the 
linear fit) at range bin i along range line j is 

/ _ -2«Sj-2j - U Sj-lj + U Sj+lj + ^Sj+2J 

\ dr s ) i j 10 Ar s 

(33) 
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The vertical divergence dw m /dz m was then com- 
puted with equation (10) or (15), depending upon the 
core criteria established in the previous section. The 
vertical divergence was then used in either the linear 
or empirical model to compute the vertical wind for 
that range bin. The computed vertical winds were 
constrained between 10 and —20 m/s to preclude 
unrealistic estimates. A flowchart of the vertical wind 
calculation process is provided in figure 16. 

The vertical wind computed from the two models 
and the true vertical wind are shown in figure 17. At 
the lower three scan altitudes, the difference between 
the two models is small, with both slightly under- 
estimating the magnitude and the spatial extent 
of the downdraft. At the upper three altitudes, 
the empirical model produces stronger downdraft 
estimates than the linear model, which continues 
to underestimate the magnitude. Neither model 
replicates the spatial extent of the downdraft at the 
upper scan altitudes. This is because of the poor 
correlation between the positive radial shear and the 
downdraft region at these altitudes. 

The error in the vertical wind estimation was 
defined as the computed value subtracted from the 
true value. The error was computed at each range 
bin of the scan for both models. The mean and the 
standard deviation of the error were then computed 
for each scan altitude. The results are shown in figure 
18 along with the error statistics that result when 
the vertical wind is neglected (assume w = 0, which 
represents the lower limit of model performance). 
The statistics confirm the qualitative assessments of 
figure 17. The mean and the standard deviation 
of the error increase with altitude. Both models 
tended to underestimate the downdraft; however, the 
empirical model did slightly better but often with a 
larger standard deviation. These results agree with 
those of reference 3, in which a similar analysis was 
conducted with the symmetrical microburst data set 
shown in figure 1. 

Vertical wind estimation with measure - 
ment error . The second part of the analysis com- 
puted the vertical winds with the simulated radar 
measurements, which included the effects of ground 
clutter and signal noise. The analysis described in 
the previous section was repeated by using the simu- 
lated radar velocity measurements shown on the right 
side of figure 12 as input to the vertical wind cal- 
culation. The mean and the standard deviation of 
the simulated radar radial velocity measurement er- 
ror were computed at each scan altitude to establish 
the fidelity of input to the vertical wind calculation. 
Figure 19 shows the mean and the standard devi- 
ation of the simulated radial velocity measurement 


error. The mean error varied from —2.0 to 2.3 with 
standard deviation from 2.6 to 6.4 m/s. 

The effect of the radar measurement error on the 
vertical wind estimation can be seen in figure 20, 
which shows the true vertical wind and the computed 
vertical wind from the two models at each scan 
altitude. When compared with figure 17, it can be 
seen that the performance of the vertical wind models 
is sensitive to the fidelity of the radar measurement. 
The magnitude of this sensitivity is presented in 
figure 21 in terms of the mean and the standard 
deviation of vertical wind estimation error which 
includes the effects of radar measurement error. The 
estimation error results with no radar measurement 
error, which were presented in figure 18, are also 
shown to establish the performance limits of the 
models. 

Vertical F-factor error. The effect of the 
wind measurement error on estimating the vertical 
component of the hazard was assessed by computing 
the vertical component of the 1-km averaged F-factor 
for each range bin. Recall that a 1-km averaged 
F-factor of 0.1 or greater is considered hazardous. 
The vertical F-factor F v was first computed for each 
range bin and then averaged over seven successive 
range bins (1050 m) along the scan line to derive the 
averaged vertical component of the hazard estimate. 
The equation for the average vertical F-factor at 
range bin i is 

.. k~i + 6 k=i + 6 

= 7 E = TV Y. < 34 > 

k—i k—i 

Figure 22 shows the true F v and that computed 
with the vertical wind estimates derived from the 
true radial winds (i.e., the vertical winds in fig. 17). 
The F v contours of figure 22 are very similar to the 
vertical wind contours of figure 17. The true F v = 
0.05 contour is highlighted in white to distinguish the 
area where the vertical contribution is at least half 
of the hazard alert threshold. Note that the area 
in which F v exceeds 0.05 can be large, even at the 
lower altitudes where the vertical wind magnitude is 
reduced. 

Figure 23 shows the F v contours computed with 
the vertical wind estimates that included measure- 
ment error effects (i.e., the vertical winds of fig. 20). 
As with figure 22, the true F v — 0.05 contour is high- 
lighted in white. The difference between the F v con- 
tours with and without measurement errors (figs. 23 
and 22, respectively) is much less than the difference 
between the vertical wind contours with and without 
measurement errors (figs. 20 and 17, respectively). 
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This difference is due to the filtering effect of 
averaging the vertical F- factor component. 

The F v error was computed in the same man- 
ner as the vertical wind error. The mean and the 
standard deviation of the F v error are shown in fig- 
ure 24 for each scan altitude. As with_the vertical 
wind errors, shown in figure 21, the^F v error in- 
creased with altitude. The mean of the F v error, with 
the measurement error effects, varied from 0.0003 to 
0.075, and the standard deviation varied from 0.015 
to 0.083. At the altitudes of primary interest for 
wind shear detection, at or below 300 m, the mean 
error was between 0.010 and 0.042 with the standard 
deviation between 0.014 and 0.028. 

The improvement in the hazard estimate achieved 
through the vertical wind estimation models was 
assessed by computing the percent improvement in 
the F v estimate. The percent improvement was 
computed relative to neglecting the vertical wind 
contribution as follows: 

Percent improvement = 100 x ^1 — 

(35) 

With the above formulation, a 100-percent im- 
provement would correspond to a perfect estimate 
from the model, and a negative percent improvement 
would indicate that the model estimate was worse 
than the neglected vertical contribution. Figure 25 
shows percent improvement in the mean error of the 
two vertical wind models with and without the mea- 
surement error effects. Figure 26 shows the corre- 
sponding percent improvement of the standard devi- 
ation. The effect of the measurement error on the 
model performance is clearly shown in these figures. 
The simulated measurement error reduces the per- 
cent improvement to about one half of that achieved 
with perfect radial velocity measurements. The mea- 
surement error effect on the standard deviation is 
much less. 

The altitude of maximum performance improve- 
ment of the models is about 300 m for both the mean 
and the standard deviation. Above 300 m, the per- 
cent improvement in the mean error gradually di- 
minishes. The percent improvement in the standard 
deviation is much less and rapidly diminishes at the 
higher altitudes. At the lower altitudes, the percent 
improvement in the mean and the standard devia- 
tion is minimized by the diminished contribution of 
the vertical wind. 

Total F-f actor error . The effect of the vertical 
wind model errors on the total hazard estimate was 
determined by repeating the F v analysis for the 


1-km F. The it term in the horizontal component 
of the total F - factor was approximated as 

U « ^ Vg (36) 

Figure 27 shows the true F and F computed with 
the vertical wind estimates derived without radar 
measurement errors. The true F — 0.1 contour is 
highlighted in white to distinguish the _areas exceed- 
ing the hazard alert threshold. The F contours of 
figure 27 are basically the same shape as the F v con- 
tours of figure 22 with about a 0.05 increase in the 
magnitude of F. This increase would indicate that 
the vertical component of the F-factor significantly 
contributes to the total F - factor for this simulated 
microburst with the aircraft at the assumed airspeed 
of 130 knots. This condition is illustrated in figure 28, 
which shows the contours of horizontal and vertical 
components of the true F-factor relative to the total. 

Figure 29 shows the true F contours and those 
computed with the vertical wind estimates that 
included measurement error effects. As with the pre- 
vious two figures, the true F = 0.1 contour is high- 
lighted in white to distinguish the hazard threshold. 
The primary effect of the measurement error on the 
F contours was an increase in the magnitude of the 
maximum and minimum values. This effect is il- 
lustrated by contrasting figure 27, which does not 
include measurement error effects, with figure 29, 
which does. This increase in the extremes is due 
to the compounding effect of the measurement error, 
which is a factor in the vertical wind estimate. The 
vertical wind estimation technique acts as an altitude 
dependent amplifier of the horizontal shear measure- 
ment. Consequently, the vertical wind estimate is 
sensitive to noise in the horizontal shear measure- 
ment, and this sensitivity increases with altitude. 

The mean and the standard deviation of the F er- 
ror are shown in figure 30 for each scan altitude. The 
mean F error was about the same as the mean F v 
error. However, because of the compounding effect of 
the measurement error, the standard deviation of the 
F error was larger. The mean of the F error, with 
the measurement error effects, varied from 0.007 to 
0.094, and the standard deviation varied from 0.029 
to 0.089. At altitudes at or below 300 m, the mean 
error was between 0.019 and 0.057, with the standard 
deviation between 0.029 and 0.049. 

Figure 31 shows percent improvement in the 
mean error of the two vertical wind models, with 
and without the measurement error effects. The 
curves for the improvement in the mean error of^F 
are similar to those for F v in figure 25. As for F ?; , 


Error mode i \ 
Error ^ =0 ) 
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the altitude of maximum performance improvement 
occurs at about 300 m. 

Figure 32 shows the corresponding percent im- 
provement of the standard deviation. The percent 
improvement in the_standard deviation of F is much 
less than that of F v . The compounding effect of 
the measurement error is clearly seen. The nega- 
tive percent improvement values at altitudes above 
200 m indicate that the standard deviation was larger 
than that from neglecting the vertical wind. The dif- 
ference between the results with and without mea- 
surement error indicates that the standard deviation 
could be improved by reducing the measurement er- 
ror below 300 m. Errors in the vertical wind models 
limit the improvement above 300 m. 

Flight Test Results 


the flowchart. The algorithm computed the residual 
of the linear least-squares fit for the radial shear. The 
residual at range bin i was defined as 



If the residual of the linear fit exceeded the residual 
threshold of 3.0 m/s, the radial shear was set to 0. 
This resulted in a vertical wind estimate of 0 for that 
range bin. 


After the radial shear and vertical wind were com- 
puted for range bins 7 through 13, the 1-km averaged 
vertical F - factor for range bin 10 was computed as 


The_flight test data analysis consisted of com- 
puting F v from the onboard Doppler radar measure- 
ment and comparing it with the in situ measurement, 
which was assumed to be the true value. The verti- 
cal wind was computed with both the linear and the 
empirical models. The results of the flight test data 
analysis were then compared with the radar simula- 
tion results. 

The onboard in situ system only measures the 
F - factor along the flight path of the airplane. There- 
fore, the radar estimates of F v can only be com- 
pared with the in situ at range bins along the flight 
path. As the airplane approaches a given point along 
the flight path, the point may have been scanned 
several times at various ranges from the airplane. 
This situation leads to a variety of possible meth- 
ods to compare the radar measurements with the 
in situ measurements. For this study, a rather sim- 
ple method was selected whereby the radar measure- 
ment at a fixed range directly in front of the airplane 
was compared with the in situ measurement of that 
airspace. The range selected was 2 km. This range 
was close enough to the airplane to minimize the ef- 
fect of microburst dynamics and airplane maneuver- 
ing between the radar and in situ measurements and 
yet provided sufficient range for the radar signal fil- 
tering algorithms to be effective. 

The radar measurements at 11 successive range 
bins were required to compute F v at a given range. 
Range bin 10 along the 0 azimuth line corresponded 
to the forward-look point at 2 km, as shown in fig- 
ure 33. The radar measurements from range bins 5 
through 15 were used in the data analysis. The radial 
shear and corresponding vertical wind estimate were 
computed with the calculation flowchart provided in 
figure 16 for range bins 7 through 13. The radar al- 
gorithm had one additional calculation not shown on 


F 


v w — 


-1 

7V 


k=13 


E w k 

k = 7 


(38) 


where V was the true airspeed of the airplane at the 
time of the radar measurement. 


To compare the in situ and the radar measure- 
ments, the time scale for the radar measurement was 
shifted by the time required for the airplane to reach 
the radar measured location. An additional 5-s shift 
was applied to correct for the lag in the in situ mea- 
surement because of its gust rejection filter (ref. 18). 
The total time shift applied to the radar measure- 
ment was 


A t = 


Unin 9.5Ar tS 


+ 5 


(39) 


For the events selected, the minimum radar range 
( r min) was 781 m and the range bin size A r s was 
144 m. The ground speed V g was that of the airplane 
at the time of the radar measurement. 


Figure 34 shows the F v of the in situ measure- 
ment and the radar estimate at range bin 10 for the 
events listed in table 4. Also shown in the figure is 
the linear correlation coefficient R of the in situ mea- 
surement and the radar estimate with the linear and 
the empirical model. The radar measurements that 
exceeded the residual threshold over the averaging in- 
terval were excluded from the correlation statistics. 
Under such conditions, radar estimate of F v defaults 
to zero. 


Figure 35 shows a summary, plotted from best to 
worst, of the correlation coefficients for the events. 
Also shown in the figure are the maximum, mini- 
mum, and average altitudes during the event. The 
best correlation between the in situ and the radar 
measurements occurred in event 548. However, this 


11 



event had a large percent of the data exceed the resid- 
ual threshold and therefore provided no estimate of 
F v . Events 438, 463, 464, 553, 555, and 573 all 
showed consistent correlation results, with values of 
R ranging from 0.778 to 0.656 with the linear model 
and 0.753 to 0.647 with the empirical model. The 
worst results were obtained from events 454, 554, 
and 556, with a negative correlation for event 454. 
These events also had the largest variation in alti- 
tude throughout the run; this indicates that perhaps 
the air mass measured by the radar was not the same 
air mass measured in situ. However, this hypothesis 
is not conclusive because event 548 yielded good re- 
sults with an altitude variation only slightly less than 
event 454. 

The estimated F v (excluding values of 0) and the 
in situ measurement for all the selected events are 
shown in figure 36. Also shown in the figure are 
the lines of perfect agreement and of the standard 
deviation of 1 and -1 about the average error. The 
average error for the linear and empirical models 
was 0.0001 and -0.0007, with standard deviations 
of 0.0087 and 0.0093, respectively. The standard 
deviation lines shown on the figure are the maximum 
of the two models. The correlation coefficient for 
the linear and the empirical model was 0.561 and 
0.557, respectively. If the events in which the altitude 
variation exceeded 200 m (events 454, 548, 554, and 
556) are excluded, then the correlation coefficient 
for the linear and the empirical model improves to 
0.659 and 0.648, respectively. The corresponding 
average error for the linear and the empirical models 
is -0.0007 and -0.0017, with a standard deviation 
of 0.0077 and 0.0083, respectively. 

Comparison of Simulation and Flight Test 

Results 

Table 4 lists the mean and the standard devia- 
tion of the F v errors from the flight test and the 
radar simulation at 300 m. The errors from the flight 
test were much less than predicted from the simula- 
tion, which could be due to a number of differences 
between the simulation and the flight test. The sim- 
ulated ground clutter environment may have been 
more severe than the clutter experienced in the flight 
test. The simulation errors were computed across the 
radar scan for a single microburst model at three dif- 
ferent times. The flight test errors were computed at 
a single range bin for a variety of microbursts events. 
The microburst model used in the simulation pro- 
duced a maximum F v value of approximately 0.20, 
which was 4 times larger than any of the flight test 
events. This larger F v resulted in the potential for 
much larger errors in the simulation. The values of 


~F V of the simulation were larger in part because the 
airspeed of the simulation (130 knots) was less than 
the flight test (210 to 230 knots). Recall from equa- 
tion (4) that the vertical F-factor is inversely pro- 
portional to the airspeed. For the same downdraft, 
an airplane traveling at 130 knots would experience 
a vertical F-factor 1.7 times greater than an airplane 
flying at 220 knots. The bottom row of table 4 lists 
the F v errors from the flight test corrected from 220 
to 130 knots. The airspeed corrected results were 
still less than the radar simulation. 

Concluding Remarks 

The objective of this study was to assess the abil- 
ity of simple vertical wind models to improve the 
hazard prediction capability of an airborne Doppler 
sensor in a realistic microburst environment. The re- 
sults indicate that, in the altitude region of interest 
(at or below 300 m), both the linear and the empirical 
vertical wind models improved the hazard estimate. 
The radar simulation study showed that the magni- 
tude of the performance improvement was altitude 
dependent. The altitude of maximum performance 
improvement occurred at about 300 m. At lower al- 
titudes, the percent improvement was minimized by 
the diminished contribution of the vertical wind. The 
performance difference between the two models was 
small. 

The results of the radar simulation study showed 
that the measurement error due to signal noise and 
clutter can significantly degrade the wind shear haz- 
ard prediction. The radar measurement errors not 
only degrade the horizontal shear hazard prediction 
but also propagate the errors to the vertical hazard 
estimate. These errors can reduce the percent im- 
provement of the vertical hazard estimate to about 
half of that achieved with perfect radial velocity 
measurements. 

The vertical hazard estimate errors from flight 
tests were less than the radar simulation results. 
This difference may be due to the lower magnitude 
of the flight test events relative to the microburst 
data sets used in the radar simulation or to an overly 
conservative simulation of the radar measurement 
error. 

The vertical hazard estimate could be signifi- 
cantly improved by reducing its sensitivity to the 
radar measurement error. This study was limited to 
processing a single radar scan to estimate the verti- 
cal wind. Methods of processing multiple radar scans 
at different elevation angles may reduce the measure- 
ment error sensitivity and improve the vertical haz- 
ard estimate. 


12 



Appendix 

Transformation Equations 

The following transformation equations relate the microburst-referenced velocities defined in a cylindrical 
coordinate system to the forward-look airborne sensor-measured velocities defined in a spherical coordinate 
system, as shown in figure 4. 

The microburst velocities (u m) v my and w m ) can be transformed to sensor-referenced velocities (u s ,v s , 
and w s ) through 


( u s \ / cos (f> s 0 sin (j> s 

Vs ] - I o 1 0 

w s / \ — sin (f> s 0 cos <j> s 

or, in matrix form, 
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Conversely, the sensor velocities can be expressed in terms of the microburst velocities as 
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The spatial derivatives of the sensor velocities, which are used in the F - factor calculation, can be derived from 
equation (A2) as 
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Equations (A4), (A5), and (A6) can be expanded by using the following geometric relationships and their 
spatial derivatives: 


/ 7 
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The partial derivatives with respect to the sensor radius r s are 
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The partial derivatives with respect to the sensor azimuth 6 S are 
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The partial derivatives with respect to the sensor elevation (p s are 
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Simplifying Assumptions 


The transformation equations are simplified by assuming the microburst is symmetrical about the vertical 


axis with no rotational velocity; therefore, 
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and 
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Spatial Velocity Gradient Equations 

Under the assumptions stated above, the spatial velocity derivatives in the microburst-centered coordinate 
system can be transformed to the sensor coordinate system as 

du s dum 2/ a a \ \ Um *2 (a a \ 2 ± 
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—— = — cos(0 m - 0 S ) cos 0 S - - — cos(0 m - 0 S ) sin 0 S 

or s or m oz rn 

du m 2 (/x a \ u m . 2(a n dw m 1 

- — cos (9 rn -9 s )-\ sin (9m-9 s )+— cos <p s sin (p s 


(A32) 


du f 

0 ^- = cos 05 1 u m sin ( - 0 S ) + Uqc sin(0oo - 6 S ) 


. (ft ft \ • 1 9 u> m , /yi \ / du m U m 

r s sin ( 1 9 m - 0 S ) sin 0 S — b cos <p s cos (0 m - 0 S ) I 

ur™ V or m r 7n , 


(A33) 


1 [Y i & u m , , du m \ 

— = cos 05 r s cos Cp s — r m tan <p s — — cos( 0 m - 0 S ) + tn m 

LV oz m dr m ) J 

+ sin <j> s r s cos <ps - r m tan <£ s - u m cos(9 m - 9 S ) - Wqo cos(0qo - #«) 

L dr m 


(A34) 


dv s . , a J du m du m 

= sin(0 m - 9 S ) r s cos <p s — r m tan 0 S — 

<705 \ dz m dr m 


(A35) 


d uj s . , f , du m du m \ 

— — = - sin 05 r s cos 0 S r m tan 0 S - — cos(0 m - 0 S ) + 

90 5 L\ 9z m 9r m J J 

+ sin 05 r s cos 0 S - r m tan 0 S - u m cos (0 m - 0 S ) - ^oo cos(0oo - 0 S ) 
L 9z m 9r m 


(A36) 


Small Angle Approximation for 0 S 

If 05 is assumed to be small (0 S < 1 ), then cos0 s « l,sin05 « 0 5 , and sin 2 0 S ^ 0; therefore, 


du s du m 2 
~drs ~ ^rm C ° S 


2 (ft a \ * /j \ - j / 

cos (0 m — 9 S ) -j sm (0 m — 9 S ) + 0 S I — 1 - — — 


/ du m dw Tn \ 

s « h )cos(^ m - 9 S ) 

V (7Z m <7r m J 


(A37) 


dvs _ /^rn _ ^m\ cos ^^ _ _ 0 s ) + s i n ( 0 m - 0 S ) (A38) 

\ OT m T xxi J OZxn 
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(A39) 
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dV m T rr, (jZjTi 


8u s 

d6 s 
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dw-n 


+ u m sin (9 m — 0 5 ) + Uoo sin (9 oo @s ) + 4> s r s ^sm(e m -0s) 
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dv s 

dO s 
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& u m -2 ta a \ . Um 2 


sm z (0 m - 0 S ) H cos z (0 m -0 s ) 

Or m r m 


dw s dw m 

= r s - — sin(0 m - e s ) - ( 
30* <9r m 


/ dUm U m \ 

^ \ r m ) 


U'm COs(0 m 0s) 'Moc COs(0(x> 0s) 
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du s du m . v , 
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(J(Ps GZm 


r m TT~“ COs(0 m - 0 5 ) - r s + Um cos(0 m - 0 S ) + Uoo COS (000 - e s 


dr Tt 


dz m 


dv s ( du m du m \ . {a a s 

-KT = )sin(0 m - 0 S ) 

c^s V dr m ) 


dw s dw r 
= r s - 
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du 

u m cos(0 m - 0 S ) - Uo o cos(0oo ~ #s) - 4>s w m + r s cos(0 m - 0 S ) + r m 

0z m 


9u r 

dr^ 


Small Angle Approximation for (f> 8 = 0 
If <f> s — 0, then 

du s du m 2/n n \ . u ^i ' 2 (a a \ 
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— = I - — ]cos(0 m - 0 S ) sm(0 m - 0 S ) 

Or 3 \Or m r m ) 
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ft7 = 6m ~ e ’ 

- r s ( - — )sin(d m - 0 s )cos(0 m - ff s ) +u m sm(d m - 0 S ) + w oo sin(0 oo - ff s ) 
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dv s 

d9 s 


= r* 
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Or m T rn 
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8w m 
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0(ps 0z m 

8v s du m . 

r s — — sm(0 m - 0 S ) 


8<j) s 8z m 
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= r x 


84> s dz- 
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(A40) 
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Table 1. Input Data for Airborne Windshear Doppler Radar Simulation Program 


Simulation parameters: 

Aircraft distance to touchdown, km Variable 0 

Aircraft velocity, knots 130 

Glide slope angle, deg 3 

Number of complete scans 6 

Time between scans, sec 3 

Roll attitude, deg 0 

Pitch attitude, deg 0 

Yaw attitude, deg 0 

Azimuth integration range/2, deg 5.0 

Azimuth integration increment, deg 0.2 

Elevator integration range/ 2, deg 5.0 

Elevator integration increment, deg 0.2 

Range integration increment, m 100 

Random number seed (0-1) 0.224 

Runway number 26 

Right (1) or left (2) 2 

Microburst and clutter: 

Microburst type (1 = 3D, 0 — 2D) 1 

Rotation of 3D microburst, deg —90 

Along track offset from touchdown, km Variable 0 

Cross track offset from touchdown, km 0.4 

Rain standard deviation, m/s 1 

Clutter standard deviation, m/s 0.5 

Clutter calculation flag (1 = ON, 0 = OFF) 1 

Discrete calculation flag (1 = ON, 0 = OFF) 1 

Reflectivity calculation threshold, dBZ —14 

Minimum reflectivity, dBZ —15 

Attenuation code (0, 1, 2) 2 

Radar parameters: 

Initial radar range, km 0.425 

Number of range cells 30 

Antenna azimuth — if no scan, deg 0 

Total number of scan lines (odd number) 15 

Azimuth scan increment, deg 3 

Antenna elevation, deg 3 

Autotilt range, km Variable 0 


Autotilt maximum elevator, deg 5 

Autotilt minimum elevator, deg — 1 

Transmitted power, W 200 

Frequency, GHz 9.3 

Pulse width, 1 

Pulse interval, fis 268.6 

Receiver noise figure, dB 4 

Receiver losses, dB 3 

Antenna type 3 

“Values shown in table 2. 
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Antenna radius, m 0.3048 

Aperture taper parameter 0.316 

Root-mean-square transmitter phase jitter, deg 0.2 

Root-mean-square transmitter frequency jitter, Hz 0 

Signal processing: 

Number of pulses 128 

Number of analog-to-digital bits 12 

Auto gain control gain factor 0.6 

Processing threshold, dB 4 

Clutter filter code (— 2 to N) 2 

Clutter filter cutoff, m/s 3 

Number of bins for F- fact or average 7 

Threshold for F-factor algorithm 0.4 

Model order autoregressive 5 

Data products: 

Alarm program data code 1 

Comma separated values plot output code 1 

Inphase quadrature phase data output code 0 

Autoregressive coefficient code 0 

Autoregressive spectra data code 0 

Fourier spectra data code 0 

Start scan number 1 

Finish scan number 6 

Start line number 7 

Finish line number 7 
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Table 2. Input Data Variables for Airborne Windshear Doppler Radar Simulation Program 


Variable 

Input data for simulation at — 

100 m 

200 m 

300 m 

400 m 

500 m 

600 m 

Aircraft distance to touchdown, km 

2.309 

4.217 

6.125 

8.033 

9.941 

11.849 

Microburst along track offset from touchdown, km 

2.892 

0.984 

-0.924 

-2.832 

-4.741 

-6.649 

Autotilt range, km 

8 

8 

0 

0 

0 

0 

Elevation angle above horizon, deg 

1.185 

0.470 

0 

0 

0 

0 


Table 3. Microburst Events Used in Flight Test Data Analysis 


Event 

number 

Location 

Radar reflectivity, 
dBZ 

Altitude, m 

In situ F 

In situ F v 

Maximum 

Average 

Minimum 

Maximum 

Minimum 

Maximum 

Minimum 

438 

Denver 

20-25 

411.2 

296.0 

243.4 

0.1167 

-0.1014 

0.0388 

-0.0192 

454 

Denver 

10-20 

559.8 

358.8 

269.1 

0.0980 

-0.1319 

0.0229 

-0.0210 

463 

Denver 

15-25 

395.0 

285.2 

248.8 

0.0881 

-0.0511 

0.0175 

-0.0091 

464 

Denver 

15-25 

453.9 

323.1 

257.0 

0.1313 

-0.0584 

0.0311 

-0.0119 

548 

Orlando 

35 45 

498.3 

328.9 

252.9 

0.1448 

-0.0807 

0.0315 

-0.0195 

553 

Orlando 

45-50 

377.5 

290.2 

235.1 

0.1470 

-0.0946 

0.0504 

-0.0094 

554 

Orlando 

30-45 

478.3 

294.3 

224.2 

0.1269 

-0.0321 

0.0252 

-0.0130 

555 

Orlando 

35-40 

439.9 

299.4 

240.2 

0.1318 

-0.0623 

0.0234 

-0.0117 

556 

Orlando 

40-45 

537.9 

340.7 

232.4 

0.0933 

-0.0687 

0.0133 

-0.0206 

573 

Orlando 

25 35 

391.9 

279.9 

222.3 

0.1202 

-0.0804 

0.0167 

-0.0134 


Table 4. Averaged F v Errors From Flight Test and Radar Simulation 
Results at 300-m Altitude 




Averaged F v errors for — 



Linear model 

Empirical model 

Test 

Mean 

Standard deviation 

Mean 

Standard deviation 

Simulation at — 

300 m, t — 49 min 

0.0188 

0.0218 

0.0100 

0.0274 

300 m, t = 50 min 

0.0313 

0.0255 

0.0196 

0.0285 

300 m, t = 51 min 

0.0415 

0.0251 

0.0326 

0.0244 

Flight 

0.0001 

0.0087 

-0.0007 

0.0093 

Flight corrected to 130 knots 

0.0002 

0.0147 

-0.0012 

0.0157 
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Figure 1. Wind vectors and F - factor contours of axisymmetric microburst simulation at four times. 
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Figure 2. F-factor contours computed without vertical winds. 
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(a) t = 49 min. 


Figure 3. Horizontal and vertical cross sections of wind vector field with radar scan area outlined. 
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Figure 3. Continued. 
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Figure 3. Concluded. 
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Figure 6. Radial variation of vertical wind of TASS generated axisymmetric microburst. (From ref. 3.) 
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Figure 9. Characteristic variation of radial shaping functions. 



Figure 10. Characteristic variation of vertical shaping functions. 
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(a) t = 49 min. 

Figure 13. Vertical wind and radial shear at each radar scan grid point. 
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Figrue 13. Concluded. 
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Figure 16. Flowchart of vertical wind calculation at range bin i along scan line j. 3 < i < Number of range 
bins —2; 1 < j < Number of scan lines. 
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(a) t = 49 min. 

Figure 17. True vertical wind and vertical wind computed from linear and empirical models with no 
measurement error at six scan altitudes for asymmetric microburst case. 
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(b) t — 50 min. 
Figure 17. Continued 





(c) t — 51 min. 
Figure 17. Concluded. 
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Scan altitude, m 

Figure 19. Mean and standard deviation of simulated radial velocity measurement error of asymmetric 
microburst case. 




(a) t = 49 min. 


Figure 20. True vertical wind and vertical wind computed from linear and empirical models with radar 
measurement errors at six scan altitudes for asymmetric microburst case. 


(b) t — 50 min. 
Figure 20. Continued 
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(a) t = 49 min. 

Figure 22. True F v and F v computed from linear and empirical models without radar measurement errors at 
six scan altitudes for asymmetric microburst case. 
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(a) t = 49 min. 

Figure 23. True F v and F v computed from linear and empirical models with radar measurement errors at six 
scan altitudes for asymmetric microburst case. 
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(b) t = 50 min. 
Figure 23. Continued. 
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Figure 26 . Improvement in F v error standard deviation 




(a) t — 49 min. 

Figure 27. The F true and F computed from linear and empirical models without radar measurement errors 
at six scan altitudes for asymmetric microburst case. 
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(a) t = 49 min. 

Figure 28. True F,F/ l> and F v at six scan altitudes for asymmetric microburst case. 




(b) t = 50 min. 
Figure 28. Continued 
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Figure 29. True F and F computed from linear and empirical models with radar measurement errors at six 
scan altitudes for asymmetric microburst case. 
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Figure 30. Mean and standard deviation of error in F for asymmetric microburst c 
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Figure 31. Improvement in F mean error. 
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(a) Event 438. 
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(b) Event 454. 

Figure 34. Estimated F and in situ measured F v time histories with corresponding R 
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Figure 35. Correlation between in situ F and radar F v for selected flight test events. 



Figure 36. Summary of in situ F v and radar F v for selected flight test events with corresponding R . 
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